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Abstract. We report observations of a new AM Herculis 
binary identified as the optical counterpart of the X-ray 
source RX J1313.2-3259, detected during the ROSAT All- 
Sky Survey (RASS). It has an orbital period of 251 min 
and is strongly modulated at optical wavelengths. The 
long-term behavior is characterized by a pronounced vari- 
ation in X-rays between the RASS and two subsequent 
pointings (decrease by a factor 40 in count rate) and by 
moderate changes in the optical brightness (up to a fac- 
tor 5). The X-ray spectrum is dominated by a soft quasi- 
blackbody component, with a smaller contribution from 
thermal bremsstrahlung. Measurements of high circular 
polarization confirm its classification as a Polar with a 
magnetic field strength of 56 MG. The average visual mag- 
nitude of RX J1313. 2-3259 is V ~ 16 m , for its distance 
we get ~ 200 pc. 
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1. Introduction 

Polars or AM Herculis binaries belong to the class of cat- 
aclysmic variables. They consist of a low mass main se- 
quence star filling its critical Roche lobe and a magnetic 
white dwarf accreting matter from its companion. The 
strong magnetic field channels the accretion flow into a 
small area on the surface of the white dwarf. In this area 
the temperature rises to values of several 100 000 K and 
shifts the emission of reprocessed radiation into the soft 
X-ray regime. Therefore an efficient way to detect these 
systems is to search the sky for (highly) variable extreme 

Send offprint requests to: hcthomas@mpa-garching.mpg.de 

* Based in part on observations collected at the European 
Southern Observatory, La Silla, Chile with the ESO/MPI 2.2m 
telescope in MPI time and with various telescopes in ESO time 
(ESO Nos. 50.6-017, 50.6-021, 54.D-0698, 55.D-0383, 56.D- 
0552, 56.D-0561, 61.D-0415). 



Fig. 1. V-band CCD image of RX J1313. 2-3259. The X- 
ray error circle from the RASS refers to the 90% con- 
fidence level. The position of the optical counterpart is 
a 2 ooo = 13 h 13 m 17 s ;l, <5 2 ooo = -32°59'12", the average vi- 
sual magnitude V ~ 16 m . 



ultraviolet or soft X-ray emitters. The great success of 
ROSAT, ASCA, and EUVE in detecting these systems 
is demonstrated by the fact, that before ROSAT only 17 
Polars were known (Cropper 1990) whilst as of today this 
number has increased to 63. Other sources of radiation in 
these systems are thermal bremsstrahlung emission from 
the accretion column, observed at higher X-ray energies 
(typically 10 to 20 keV), cyclotron emission from the ac- 
cretion column due to the strong magnetic field, often 
dominating the optical regime, and radiation from the sec- 
ondary star, strongest in the infrared and sometimes only 
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Table 1. List of observations. Columns denote: (1) and (2) the time of the observation, (3) and (4) the telescope 
and instrument used, (5) the type of observation (spec: spectrophotometry, point.: X-ray pointing, phot.: photometry, 
sp.pol.: spectropolarimetry), (6) the spectral range or the filters used, (7) the FWHM resolution (spectroscopy only), 
(8) the number of spectra or V filter data obtained, (9) the time span for the observation, and (10) the exposure time 
for the individual measurements. 
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being detectable, if the system is in a state of low accre- 
tion. Detailed discussions of these systems can be found in 
the books by Warner (1995) and Campbell (1997). Early 
reviews on the basis of the new ROSAT data have been 
presented among others by Beuermann & Thomas (1993), 
Watson (1994), Beuermann & Burwitz (1995), and Schwo- 
pe (1995), a more recent one by Beuermann (1997). The 
distribution of Polars in the solar neighborhood was inves- 
tigated by Thomas & Beuermann (1997). An identification 
program based on a complete sample of the brightest soft 
X-ray sources from the RASS (Voges et al. 1999) at high 
galactic latitudes (Thomas et al. 1998) led to the detection 
of RX J1313.2-3259. 



2. Observations 

Starting 12 January 1991 the position of RX J1313. 2-3259 
in the sky was scanned by the ROSAT XRT with the 
PSPC as detector. During 20 satellite orbits with a total 
exposure time of 485 s the object was detected with a 
mean count rate of 1.8 cts s -1 and a hardness ratio HRI = 



-0.85Q. On a finding chart produced from the COSMOS 
scans of the SERC-J plates at ROE (Yentis et al. 1992) a V 
= 16 m star was found to be the likely optical counterpart, 
at a distance of ~ 3" from the X-ray position (Fig. 1). 
A low resolution spectrum taken on 31 August 1991 at 
the ESO/MPI 2.2m telescope showed strong Balmer and 
helium line emission thus confirming this identification. A 
list of all observations collected until now can be found in 
Table 1. 

2.1. X-ray photometry and spectroscopy 

From the RASS data a photon event table was extracted 
which covered 50 x 50 arcmin 2 around the X-ray position 
of the source. Using the EXSAS software package provided 
by the MPE Garching (Zimmermann et al. 1994) for the 
extraction of source photons in a circle of 250" radius cen- 
tered on the source and background photons in a circle of 
400" radius in scan direction we obtained the light curve 

1 HRI = (H - S)/(H + S) with H and S the count rates 
in the hard and soft energy intervals 0.5 — 2.4 keV and 0.1 — 
0.4 keV, respectively. 
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shown in Fig. 2 and the mean spectrum inserted in Fig. 4. 
Phasing of the light curve has been obtained from Eq. (2) 
based on optical observations (see Sect. 3.2), so phase zero 
corresponds to the inferior conjunction of the secondary 
star. The error in the period results in a phase error be- 
low 0.01, therefore the cycle count is correct. The data are 
sampled from nine different orbital cycles of the binary. 
The mean count rate obtained from the light curve is 1.94 
ctss -1 . In one ROSAT orbit a count rate far above the 
average (6.5 ctss -1 ) was measured. We consider this as a 
singular event not typical for the orbital variation. With- 
out this point the mean count rate drops to 1.68 ctss -1 . 
In order to show that the light curve is not dominated 
by changes between binary orbits we have used two dif- 
ferent symbols in Fig. 2 for the first and second half of 
the observation. A period search without the point at 6.5 
ctss" 1 revealed a likely period of 254 ± 7 min. It should 
be noted that at no phase did the count rate drop to zero, 
so the accretion area giving rise to the X-ray emission 
never completely vanishes from view behind the horizon 
of the white dwarf. The alternative of a second accretion 
area contributing to the X-ray flux seems to be ruled out 
by the results of our spectropolarimetry, which shows no 
change of sign in the circular flux (see Sect. 2.3). 

In two subsequent pointings with the ROSAT PSPC 
(July 1992) and HRI (July 1994) as detectors the source 
was observed with count rates of 0.036 PSPC ctss -1 
(hardness ratio HRI = 0.28) and 0.010 HRI ctss -1 (cor- 
responding to ~ 0.06 PSPC ctss -1 ). This is a reduction 
in the count rate by factors 50 (1992) and 30 (1994) com- 
pared to the RASS. Phase binning of these data resulted 
in the lightcurves shown in Fig. 2. Again we note that the 
site of the X-ray emission always remains in view of the 
observer. 

Fitting the RASS spectrum with a blackbody and a 
thermal bremsstrahlung component together with possi- 
ble interstellar absorption resulted in a blackbody tem- 
perature of 58 eV and a column density of 9 • 10 19 
H- atoms cm -2 . The uncertainties for these parameters are 
depicted in Fig. 3. The (unabsorbed) contributions of the 
two components to the total flux in the ROSAT window 
(0.1 to 2.4 keV) amount to 1.14 • 10 -11 crgcm -2 s -1 and 



0.07- 10" 11 erg cm" 



respectively. For the temperature 



of the thermal bremsstrahlung component we assumed a 
value of 10 keV, since this cannot be derived from the spec- 
tral data. The resulting unabsorbed spectra are shown in 
Fig. 4 as dotted lines. Separating the count rate into the 
two components one obtains 1.70 ctss -1 for the blackbody 
component and 0.07 ctss -1 for the thermal bremsstrah- 
lung component. The flux ratio Fbrems/^bb in the ROSAT 
band is 0.06, integrated over all frequencies it increases to 
0.16. The uncertainty in these values is about a factor 2. 

The spectral fit to the data from the PSPC pointing 
gave a blackbody temperature of 50 eV, an unabsorbed 
flux of 4.7 • 10 -14 crgcm -2 s -1 for the blackbody compo- 
nent and 9.1 • 10 -13 ergcm -2 s -1 for the thermal brems- 




phase 

Fig. 2. X-ray light curve of RX J1313. 2-3259. The upper 
panel shows the source during the RASS (first half of ob- 
servation: filled circles, second half: open circles) together 
with the measured background (crosses), the lower panels 
during two subsequent pointings with the PSPC and the 
HRI as detectors, respectively. For the preliminary period 
determination the single high datum at 6.5 ctss -1 has 
been ignored. The phases are computed from Eq. (2) in 
Sect. 3.2. 



strahlung component, with large uncertainties and a high 
correlation between flux and temperature of the black- 
body fit (see Sect. 3.4). Here we again fixed the tempera- 
ture for the thermal bremsstrahlung component to 10 keV, 
and the column density to the value obtained during the 
RASS. The data are also shown in Fig. 4. 

2.2. Optical photometry 

Sequences of direct images in different filters were taken 
at several epochs (see Table 1). To demonstrate the long- 
term variability of RX J1313. 2-3259 we have plotted in 
Fig. 5 the mean values of V in the different runs together 
with the range of variability defined by the magnitudes 
which bracket 90% of the measurements in the correspond- 
ing observing run. Added to this plot are V-magnitudes 
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Fig. 3. Confidence levels for the X-ray spectral fit of 
RX J1313.2-3259. The 1, 2, and 3cr ranges for the fit pa- 
rameters blackbody temperature and column density are 
shown together with the fit result (cross). 
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Fig. 4. Overall observed spectrum of RX J1313.2-3259. 
The optical spectra shown at low frequencies were ob- 
tained in Feb. 93 (high state) and in Dec. 95 (low state). In 
the X-ray regime the RASS data are shown as circles, the 
PSPC pointing data as crosses. The dotted lines give the 
unabsorbed blackbody and bremsstrahlung components of 
the X-ray spectral fit to the RASS data. 



from spectroscopy, obtained by folding the mean spectrum 
of each observing run with the sensitivity of the V-filter 
(a flux of 3.64 • I0~ 9 ergcm _2 s _1 A _1 corresponds to zero 
magnitude, see Bessell 1979). The dotted lines indicate the 
times of the X-ray observations (RASS, PSPC pointing, 
and HRI pointing). Typically RX J1313.2-3259 is found 
at V ~ 16 m and brightened only for short episodes near 
HJD 2 449 036 (Feb. 93) and HJD 2 450 936 (May 98). 

The orbital variation of RX J1313.2-3259 in the V- 
band is displayed in Fig. 6. Different symbols identify the 
different observing runs which are listed at the bottom of 
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Fig. 5. Long-term light curve of RX J1313.2-3259. Photo- 
metric observations are marked with filled dots (mean val- 
ues) and a vertical bar (90% range), magnitudes derived 
from mean spectra (less accurate) by stars. The dotted 
lines show the times of X-ray observations (RASS, PSPC, 
HRI). 



the figure. The lightcurves are clearly of different character 
in the bright (Feb. 93, May 98) and faint (all other) states 
of RX J1313.2-3259. While the bright states show only 
one minimum per orbit, the variation in the faintest state 
(Feb. 95) almost follows a sinusoidal variation with two 
minima per orbit. With increasing brightness the flux first 
increases outside the primary minimum (near phase zero) 
and starts to deviate from its value around the primary 
minimum for the bright states only. 

It is worth mentioning that the observations in Feb. 93 
were collected in six subsequent nights with two different 
instruments and that during this high state episode the 
orbital-averaged flux decreased by 0.07 mag per day, sug- 
gesting a characteristic duration of the high state of a 
few weeks. The orbital light curve for the Feb. 93 data in 
Fig. 6 (diamonds) was obtained by subtracting a linear 
trend from the data. The fitted light curve (dotted line) 
represents a sinusoidal fit to the detrended data. 

2.3. Optical spectroscopy and spectropolarimetry 

Time sequences of spectra at different resolutions were ob- 
tained at 9 epochs as listed in Table 1 and displayed in 
Fig. 5 (stars). We first display in Fig. 7 the low resolution 
spectra from Feb. 93 (brightest state) and from July 95 
(faintest state) . Both spectra clearly show the contribution 
of a late- type star to the flux at the red wavelengths, iden- 
tifiable through its TiO absorption troughs. The emission 
lines of hydrogen and helium are visible in both spectra, 
with He ii A4686 strongly reduced in the fainter spectrum. 

Complete coverage of the orbit in one night could be 
achieved for the observation run in March 97 (interme- 
diate state) only. Medium resolution spectra in the blue 
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Fig. 6. Phase-folded lightcurves at different observing times as given at the bottom of the plot. For the dotted line 
see text in Sect. 2.2. The phases are computed from Eq. (2) in Sect. 3.2. 



Table 2. Results from double Gaussian fits. The columns 
list the fitted emission line, the velocity amplitudes in 
km/s and phases for blue-to-red zero crossing both for the 
narrow and the broad component and the line flux ratio 
of minimum to maximum for the narrow component. 
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498 
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HeiiA4686 
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0.26 


Can K 
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0.972 
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(3600-5200 A) and in the red (6440-8360 A) were an- 
alyzed for Doppler shifts in their emission and absorp- 
tion lines. The most accurate results were obtained for 
the red spectra. Fitting the Ha emission line with a dou- 
ble Gaussian profile splits the line into a narrow (average 
FWHM 5.9 A, unresolved) and a broader (average FWHM 



17 A) component. The resulting radial velocities are plot- 
ted in Fig. 8, upper panel. The average radial velocities 
of three absorption lines (Nai A8183, 8195 and KiA7699) 
are shown in the same panel. In the lower panels of Fig. 8, 
the fluxes in the narrow emission component of He i A6678 
and the absorption line of K i A7699 are displayed together 
with results from irradiation computations (see Sect. 3.2). 
The measured fluxes of the absorption line are very sen- 
sitive to the assumed level of the continuum near this 
line and may be in error by up to 30%. We also mea- 
sured the fluxes in the narrow emission components of 
Ha, H/3, H 7 , m, HenA4686, HeiA4026, and Can K by 
fitting double Gaussians to the profiles. In Table 2 the 
velocity amplitudes and phases for the blue-to-red zero 
crossings are given for the two components together with 
the line flux ratios (minimum to maximum) of the narrow 
components. The phases were derived from the ephemeris 
given in Sect. 3.2. The minimum emission line flux stays 
finite for all of these lines. In this respect, all lines behave 
similarly to HeiA6678 in Fig. 8. 

During the observation in Dec. 95 the source was in 
a low state, which allows for a separation of the spectral 
flux into the different contributions from the M-star, the 
accretion area, and the white dwarf. The results of this 



6 



H.-C. Thomas et al.: RX J1313. 2-3259, a long-period Polar 




4000 



5000 



6000 7000 
Wavelength [A] 



8000 



9000 



Fig. 7. Spectra at low resolution from two different ob- 
serving runs. The brighter spectrum was taken in Feb. 93, 
the other is the mean of six spectra obtained in July 95. 

analysis will be discussed in Sect. 3.3. The spectra taken 
during another low state in Jan. 92 are not suited for this 
kind of analysis because they cover less than half an orbit 
and their flux calibrations are unreliable. 

For the observation in May 98 (high state) EFOSC 2 
was equipped with the Wollaston prism and a quarterwave 
plate, using Grism B300. The sequence of 2 x 9 spectra 
covers only 60% of the orbital period. Circularly polar- 
ized fluxes have been obtained by taking the difference 
between the two spectra produced by the Wollaston prism 
(Fig. 9). All spectra show circular polarization of negative 
sign only. We find, that two maxima of the polarized flux, 
at ~ 5000 A and ~ 6600 A, are present in all spectra ex- 
cept near phase where the viewing angle is smallest with 
respect to the axis of the accretion funnel. In addition, a 
third maximum of the polarized flux occurs at ~ 3950 A 
over the restricted phase interval of 0.3 to 0.5. We will 
argue below that these are the 3 rd , 4 th , and 5 th harmonic 
in a field of about 56 MG and that the different phase 
behavior arises from optical depth and geometric effects 
(Sect. 3.3). 

Spectra taken during the other observing runs have 
been analyzed for radial velocities, flux contribution from 
the secondary star, and cyclotron emission. Because of 
lower resolution, incomplete orbital coverage, or less favor- 
able observing conditions they mainly helped to reinforce 
and confirm the results obtained from the March 97 data 
and are important for excluding possible alias periods. 

3. Results 

3.1. A distance estimate 

In all low resolution spectra the TiO absorption bands 
from the secondary are easily recognized. We therefore use 
the calibration method dccsribed in Beucrmann & Welch- 
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Fig. 8. Line properties extracted from spectra taken in 
Mar. 97. Upper panel: radial velocities averaged from mea- 
surements of three absorption lines (stars), from the nar- 
row component of the Ha emission line (filled circles), 
and the broad component of the Ha emission line (open 
circles). The solid, dashed, and dotted lines give the si- 
nusoidal fits to the data. Middle panel: normalized line 
flux of the He i A6678 narrow emission component together 
with results from irradiation computations for three dif- 
ferent inclinations (solid: 66°, dashed: 52°, dotted: 31°, 
see Sect. 3.2). Lower panel: normalized line flux of the 
KiA7699 absorption line together with results from irra- 
diation computations for the same inclinations. 



hold (1999) to determine the surface flux from the flux 
difference between 7165 and 7500 A. In July and Dec. 95 
RX J1313. 2-3259 was in a low state, so we used the flux 
differences as observed. For the Mar. 97 data we first es- 
timated the gradient of the cyclotron emission between 
7165 A and 7500 A by subtracting suitably scaled spectra 
of the M-dwarf Gl 207.1 and corrected the measured flux 
differences for that. We also tried to fit spectra of other 
M-dwarfs (e.g. Gl 205, Gl 352) but with less satisfactory 
results. From these 35 spectra we obtained a mean value 
of 9.7 ± l.l-10 _16 ergcm~ 2 s _1 A _1 for the flux difference. 
This value is reduced by 2% if we exclude all spectra with 
phases between 0.25 and 0.75, i.e. looking at the unillumi- 
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Fig. 9. Circularly polarized flux extracted from spectra 
taken in May 98. The grayscale plot displays fluxes be- 
tween (white) and — 6T0~ 16 crgcm _2 s _1 A _1 (black), 
orbital phases are given on the right-hand side. The image 
is smoothed in order to remove the contributions from the 
emission lines. 

nated backside of the secondary. In the following we use a 
value of 9.5-I0 _16 ergcm _2 s _1 A _1 . 

Filling its critical Roche lobe the secondary must have 
a mean density of p/pe — 4.46 for a system period of 
4.19 h. This value depends slightly on the mass ratio but 
the distance changes by no more than 2% for white dwarf 
masses between 0.5 and 1.2 M©. Assuming the secondary 
to be a main-sequence star the stellar models of Baraffc 
et al. (1998) suggest a mass of 0.51 M Q and a radius of 
0.486 R Q , with a slight dependence on mass ratio. How- 
ever, the secondary, due to its mass-loss history, may be 
out of thermal equilibrium and therefore have a lower 
mean density than a main-sequence star of the same mass. 
Comparing the low-state spectra to those of several M- 
stars we estimate the spectral type to be M2.5 ±0.5. Since 
the spectral type does not change drastically for stars out 
of thermal equilibrium (Kolb & Baraffe 1999), we obtain 
an estimate for the mass of the secondary from a compar- 
ison with models of Baraffe et al. (1998): M 2 /M© ~ 0.45. 
Using the relation between spectral type and the surface 
brightness Fxio by Beuermann & Weichhold (1999) we 
arrive at a distance of 

200±17pc- (M 2 /0.45M Q ) 1/3 . (1) 
3.2. System parameters 

The orbital period of RX J1313. 2-3259 was obtained from 
a sinusoidal fit to the radial velocity measurements of the 
absorption lines Nai A8183,8195 and KiA7699, obtained 
in Aug. 92, Feb. 93, and Mar. 97. Possible aliases have 
been checked against other radial velocity measurements 
obtained at many different epochs (see Table 1). The zero 
point for the resulting ephemeris is defined as the blue-to- 
red zero-crossing of the radial velocity curves. The fits to 
the narrow emission line components of the Mar. 97 data 
give a slightly earlier zero point than the absorption line 



data (difference in phase: 0.06 ± 0.03) and we have used 
the mean of both data sets. The resulting ephemeris is 
(errors of the last digits are given in brackets): 

T = HJD 2 448 632.4435(43) + 0.174 592 09(15) E (2) 

The relatively large error in the zero point reflects the dif- 
ference in determining the bluc-to-red zero-crossing from 
both the narrow emission and the absorption line mea- 
surements. The total time span between our first and last 
observation amounts to 2667 days or 15276 cycles. Thus 
the maximum error in the phasing calculated from the er- 
ror for the period is 0.013, which excludes any errors in 
the cycle count. 

The measurements shown in the two lower panels of 
Fig. 8 clearly demonstrate that the narrow emission line 
flux is stronger and the absorption line flux weaker than 
average when the illuminated side of the secondary is 
in view of the observer. Such a behavior can be under- 
stood with the irradiation model of Beuermann & Thomas 
(1990). Using the amplitude of 185 ±9 km/s measured for 
the absorption lines the model determines the inclination 
angle as a function of mass ratio. We have assumed that 
the contribution to the absorption line flux from the illu- 
minated side is negligible. For primary masses of 0.5 and 
1.1 M© we obtain inclinations of 52° and 31°, respectively, 
taking 0.45 M© for the mass of the secondary. The corre- 
sponding velocities for the emission lines are then 51 and 
98 km/s. Looking at the measured values (Table 2) this 
would favor the low inclination and so a high primary 
mass. A low inclination also roughly fits the observed line 
fluxes, an example is shown in Fig. 8 (middle panel, dotted 
line). For the inclination of 66° derived below (Sect. 3.3) 
this is not the case. But if we assume that a substan- 
tial contribution to the line flux comes from the unillu- 
minated side of the secondary the situation changes. To 
reproduce the observed average amplitude of 102.6 km/s 
for the narrow emission lines, ~ 55% of the line flux (aver- 
aged over the surface of the secondary) must be generated 
by illumination. We then obtain the line flux variation dis- 
played in Fig. 8 (middle panel, solid line), which fits the 
observed variation best, except for those measurements 
where the radial velocity curves of broad and narrow com- 
ponent cross each other and therefore a separation into 
the two components is difficult to achieve. While this re- 
sult depends only weakly on the assumed inclination it 
shows that the observed line flux and velocity amplitude 
for the narrow emission are consistent with a high incli- 
nation if part of the emission line flux is generated on the 
unilluminated side of the secondary. 

3.3. Cyclotron emission 

For the data taken during a low state in Dec. 95, we rep- 
resent the spectral flux in the 14 individual low-resolution 
spectra by the sum of a Rayleigh-Jeans spectrum, an M- 
star spectrum, and a cyclotron spectrum. This was done 
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Fig. 10. Cyclotron flux extracted from spectra taken in 
Dec 95. The model fits are discussed in Sect. 3.3, orbital 
phases are given on the right-hand side. 



by adjusting the Rayleigh- Jeans to the flux between 4040 
and 4080 A and the M-star template to the flux differ- 
ence between 7165 and 7500 A. (As discussed in Sect. 3.1 
the best M-star template is that of the M3 star Gl 207.1). 
Subtracting these two contributions we are left with the 
cyclotron spectra shown in Fig. 10. The variation of the 
Rayleigh- Jeans component in the V-band is roughly sinu- 
soidal and can be fitted with a mean of 3.1±0.1 and an am- 
plitude of 0.3±0.1, both in units of 10~ 16 erg cm _2 s _1 A _1 , 
the maximum occuring at phase 0.01 ± 0.4. Such a vari- 
ation indicates heating of an area around the accretion 
spot (see Sect. 3.3). The flux of the M-star displays max- 
ima at phases near quadrature. The origin of the these 
variations must be due to ellipsoidal modulation from the 
M-star. The flux ratio between phases 0.5 (primary mini- 
mum) and 0.25/0.75 (maxima) amounts to 0.85 ± 0.02. 

To model the Dec. 95 variation of cyclotron flux with 
orbital phase we adapted the emission of an isothermal 
homogeneous plasma slab (Barrett & Chanmugam 1985) 



to the data. The input parameters were the magnetic field 
strength B, the plasma temperature T, the inclination i 
of the system, the angle (3f between the rotation axis and 
the field direction, the angles f3 s and ip s describing the 
location of the accretion region (relative to the rotation 
axis and the line connecting the two stars), and the di- 
mensionless thickness A and the area A s of the accretion 
region. From the phasing of the broad emission lines we 
took the azimuthal angle between field direction and di- 
rection to the companion rpf to be 17°. The fit shown in 
Fig. 10 required B = 56 MG, T = lOkeV, i = 66°, /?/ = 
7°, a = 24°, tp s = 18°, A = 26, and A s = 8.2 • 10 16 cm 2 
(for a distance of 200 pc). It was obtained in a iterative 
procedure starting with a coarse grid of cyclotron spectra 
which allowed to fix B, T, i, and /?/ and then refining the 
grid for the other parameters. The poor fit at phase 0.915 
may be the result of absorption by the accretion stream 
which is in front of the accretion area around phase 0.95. 
The rather high inclination can only be made consistent 
with the results from the irradiation model for a white 
dwarf mass of 0.39 M Q (see Sect. 3.5). The values of /?/ 
and P s are in conflict with the assumption of a pure dipo- 
lar field configuration because that requires (3f > (3 S . Also 
this simple model does not reproduce the spectral shape 
at short wavelengths, probably due to a more complicated 
accretion geometry than used here. 

Using the parameters derived above we computed the 
gravity darkening in the Roche geometry. For the ratio of 
primary minimum to maximum we obtained a value be- 
tween 0.87 (without limb-darkening) and 0.83 (50% limb- 
darkening), in agreement with the value for the ellipsoidal 
variation deduced above. 

For the values of i, /3 S , and i[) s the sinusoidal varia- 
tion of the Rayleigh- Jeans component can be explained 
assuming a constant contribution from the white dwarf 
photosphere plus a varying contribution from the accre- 
tion area. The fluxes in the V-band then amount to 2.8 and 
0.85, respectively, both in units of 10~ 16 ergcm _2 s _1 A _1 . 
At the derived distance the constant flux gives an abso- 
lute magnitude of My = 11. 3 m , corresponding to a pho- 
tospheric temperature of ~ 14 000 K. The flux ratio of the 
variable to constant component in an accretion area occu- 
pying a fraction / on the white dwarfs surface requires a 
temperature which is a factor 0.74 /~ 25 higher than the 
temperature of the white dwarf, in reasonable aggreement 
with the results of Gansicke et al. (1999). 

We also tried to fit the spectropolarimetric data from 
May 98, which do not suffer from the uncertainties intro- 
duced by the subtraction of other flux contributions. But 
because of the incomplete orbital coverage we can only 
state that magnetic field strength, inclination, and field 
direction are similar, while the thickness A must be some- 
what larger to produce the high circularly polarized flux 
at 3950 A around phase 0.4 (see Fig. 9). We plan to re- 
peat these observations, which will then allow us to check 
the system parameters derived above. Cyclotron spectra 
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extracted from the Mar. 97 data show similar variations 
with phase as the data of Dec. 95, but the removal of other 
flux contributions introduced too large an uncertainty, so 
we did not try to fit cyclotron spectra to that dataset. 

These considerations also allow us to qualitatively un- 
derstand the behavior of the V-flux variations (Fig. 6). 
In a low state both the cyclotron flux and the ellipsoidal 
variations of the secondary cause a light curve with two 
maxima during one orbit, while in a high state the con- 
tribution of the secondary is negligible and the cyclotron 
flux contribution most likely comes from a more extended 
accretion area. The minimum around phase 0.5 which re- 
sulted from a viewing angle close to 90° will then be filled 
in because of contributions from other parts of the accre- 
tion area at lower viewing angles. 

3.4- X-ray emission 

We now turn to the X-ray light curves and spectra ob- 
tained during the RASS and the two ROSAT PSPC and 
HRI pointed observations. Since the inclination is larger 
than the field line direction it is expected that at least 
part of the X-ray flux is absorbed by the stream of mat- 
ter towards the white dwarf. This explains the low flux 
around phase 0.95 both in the RASS and to some degree 
in the HRI pointing. The reduction in the count rate by a 
factor 50 between the RASS and the PSPC pointing and 
the corresponding change in the hardness ratio HRI from 
—0.96 to 0.28 could partially be the result of a lower black- 
body temperature. Although the best fit temperature for 
the pointing is 50eV (x 2 /d.o.f. 22/10), with a reduction 
of the blackbody flux by a factor 230 compared to the 
RASS, the fit allows for much higher blackbody fluxes if 
the temperature would be lower: at 15 eV, the total black- 
body flux decreases by a factor 14 only compared to the 
RASS (x 2 / d -°- f - 27/11). This is of the same order as the 
reduction in cyclotron flux between the observations in 
Feb. 93 and Dec. 95 (factor 7.4). The contribution of the 
bremsstrahlung flux to the PSPC pointed data is reduced 
by a factor 2.3 compared to the RASS. 

3.5. Masses of the components 

The masses determined above are in conflict with present 
understanding of stable mass transfer. For donor stars 
with masses below ~ 0.5 M mass transfer is dynamically 
unstable for mass ratios M 2 /Mi > 0.7 (Webbink 1985). 
The lower limit of 0.38 M Q for the mass of the secondary 
(Sect. 3.1) also does not fulfil the criterion for stability. 
Even if we assume an evolved secondary as in Kolb & 
Baraffe (1999) we arrive at masses of 0.31 M Q for the white 
dwarf and 0.25 M Q for the M star, again violating the sta- 
bility criterion. The mass ratio is determined from the 
observed absorption line velocity and the inclination. So 
we turn the question around and ask, for which inclination 
would the stability criterion be satisfied? Taking the de- 



rived mass of the secondary we then need a primary mass 
of0.64M Q and an inclination of 45.3°. That implies a vari- 
ation of the viewing angle around this value, which in turn 
causes cyclotron emission lines to be shifted to the blue as 
compared to our data in Fig. 10. Especially above 8000 A 
the spectral flux should strongly increase with wavelength 
for all phases due to the presence of the next lower har- 
monic, which in Fig. 10 only marginally shows up around 
phase zero. Therefore with our data we can see no remedy 
to the present situation. So the stability problem makes 
RX J1313. 2-3259 an interesting system for future observa- 
tions in four aspects: a) the contribution of the secondary 
star to the spectrum, best observable during a low state, 
b) the absorption line flux variation relevant for modelling 
the illumination, c) the degree of circular polarization as 
a function of orbital phase, and d) the white dwarf contri- 
bution to the spectrum, best observable in the ultraviolet 
regime. 

Finally, we note that constraining the mass of the 
white dwarf is important for our understanding of the for- 
mation of cataclysmic variables, because with the period 
determined above models of common envelope evolution 
predict that the mass of the white dwarf should exceed 
0.6 M (de Kool 1992). 

3.6. Mass accretion rates 

Table 3. Integrated fluxes from the accretion area and 
mass accretion rates. The units are 10 -11 ergcm~ 2 s _1 for 
the fluxes and 10 -11 M yr -1 for the mass accretion rates. 
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4.3 


- 7.7 
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For an estimate of the mass accretion rate one has to 
sum up all flux contributions from the accretion area. Al- 
though the optical and X-ray observations are not simul- 
taneous we take the optical observations in Feb. 93 and 
the X-ray observations during the RASS as representative 
of the high state and those of Dec. 95 in the optical and 
July 92 in X-rays of the low state. To extend the obser- 
vations to the whole frequency range we used the model 
fits for cyclotron emission, blackbody radiation, and ther- 
mal bremsstrahlung, respectively. The main uncertainty 
in these values is caused by the uncertainty in the black- 
body temperature. Therefore we took the la upper and 
lower limits of 68 and 41 cV (see Fig. 3) for the RASS data 
and assumed temperature limits of 15 and 50 eV for the 
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PSPC data (see Sect. 3.4), to provide a range for the de- 
rived fluxes. The results are summarized in Table 3. From 
the total fluxes we computed the mass accretion rates for 
a distance of 200 pc, a white dwarf mass of 0.4 M , and 
a white dwarf radius of 1.08 • 10 9 cm. 

4. Conclusions 

It has been shown that RX J1313. 2-3259 belongs to the 
class of magnetic cataclysmic variables called AM Hcr- 
culis binaries or Polars. It has a period of 251.4 min, 
which places it above the period gap, and a distance of 
~ 200 pc. From the occurrence of cyclotron humps a mag- 
netic field strength of 56 MG for the main accreting pole 
has been deduced. The inclination is i ~ 66° implying 
a white dwarf mass of ~ 0.4 M Q , with a likely mass for 
the secondary of 0.45 M Q , corresponding to a mass ratio 
q ~ 1.1. A rough estimate of the mass accretion rate leads 
to 6-10" 11 M yr" 1 for the high state and 6-10" 12 M yr" 1 
for the low state. The rather low mass of the white dwarf 
required for reproducing the observed cyclotron emission 
is in conflict with our present understanding of stable mass 
transfer and the evolutionary history of cataclysmic vari- 
ables, so further observations of this system are of great 
interest. 
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